Extensive androgen receptor enhancer heterogeneity in primary prostate cancers underlies transcriptional diversity and metastatic potential.

Androgen receptor (AR) drives prostate cancer (PCa) development and progression. AR chromatin binding profiles are highly plastic and form recurrent programmatic changes that differentiate disease stages, subtypes and patient outcomes. While prior studies focused on concordance between patient subgroups, inter-tumor heterogeneity of AR enhancer selectivity remains unexplored. Here we report high levels of AR chromatin binding heterogeneity in human primary prostate tumors, that overlap with heterogeneity observed in healthy prostate epithelium. Such heterogeneity has functional consequences, as somatic mutations converge on commonly-shared AR sites in primary over metastatic tissues. In contrast, less-frequently shared AR sites associate strongly with AR-driven gene expression, while such heterogeneous AR enhancer usage also distinguishes patients' outcome. These findings indicate that epigenetic heterogeneity in primary disease is directly informative for risk of biochemical relapse. Cumulatively, our results illustrate a high level of AR enhancer heterogeneity in primary PCa driving differential expression and clinical impact.

Optimizing Efficiency and Motility of a Polyvalent Molecular Motor.

Molecular motors play a vital role in the transport of material within the cell. A family of motors of growing interest are burnt bridge ratchets (BBRs). BBRs rectify spatial fluctuations into directed motion by creating and destroying motor-substrate bonds. It has been shown that the motility of a BBR can be optimized as a function of the system parameters. However, the amount of energy input required to generate such motion and the resulting efficiency has been less well characterized. Here, using a deterministic model, we calculate the efficiency of a particular type of BBR, namely a polyvalent hub interacting with a surface of substrate. We find that there is an optimal burn rate and substrate concentration that leads to optimal efficiency. Additionally, the substrate turnover rate has important implications on motor efficiency. We also consider the effects of force-dependent unbinding on the efficiency and find that under certain conditions the motor works more efficiently when bond breaking is included. Our results provide guidance for how to optimize the efficiency of BBRs.

Functional mapping of androgen receptor enhancer activity.

BACKGROUND: Androgen receptor (AR) is critical to the initiation, growth, and progression of prostate cancer. Once activated, the AR binds to cis-regulatory enhancer elements on DNA that drive gene expression. Yet, there are 10-100× more binding sites than differentially expressed genes. It is unclear how or if these excess binding sites impact gene transcription. RESULTS: To characterize the regulatory logic of AR-mediated transcription, we generated a locus-specific map of enhancer activity by functionally testing all common clinical AR binding sites with Self-Transcribing Active Regulatory Regions sequencing (STARRseq). Only 7% of AR binding sites displayed androgen-dependent enhancer activity. Instead, the vast majority of AR binding sites were either inactive or constitutively active enhancers. These annotations strongly correlated with enhancer-associated features of both in vitro cell lines and clinical prostate cancer samples. Evaluating the effect of each enhancer class on transcription, we found that AR-regulated enhancers frequently interact with promoters and form central chromosomal loops that are required for transcription. Somatic mutations of these critical AR-regulated enhancers often impact enhancer activity. CONCLUSIONS: Using a functional map of AR enhancer activity, we demonstrated that AR-regulated enhancers act as a regulatory hub that increases interactions with other AR binding sites and gene promoters.

Substrate stiffness tunes the dynamics of polyvalent rolling motors.

Nature has evolved many mechanisms for achieving directed motion on the subcellular level. The burnt-bridges ratchet (BBR) is one mechanism used to achieve superdiffusive molecular motion over long distances through the successive cleavage of surface-bound energy-rich substrate sites. This mechanism has been associated with both nanoscale and microscale movement, with the latter accomplished through polyvalent interactions between a large hub (e.g. influenza virus) and substrate (e.g. cell surface receptors). Experimental successes in achieving superdiffusive motion by synthetic polyvalent BBRs have raised questions about the dynamics of their motility, including whether rolling or translation is better able to direct motion of microscale spherical hubs. Here we simulate the three-dimensional dynamics of a polyvalent sphere moving on and cleaving an elastic substrate. We find that substrate stiffness plays an important role in controlling both the motor's mode of motility and its directional persistence. As we tune lateral substrate stiffness from soft to stiff we find there exists an intermediate value that optimizes rolling behaviour. We also find that there is an optimal substrate stiffness for maximizing persistence length, while stiffness does not influence as strongly the superdiffusive dynamics of the particle. Lastly, we examine the effect of substrate density, and show that softer landscapes are better able to buffer against decreases in substrate occupancy, with the spherical motor maintaining superdiffusive motion more on softer landscapes than on stiff landscapes as occupancy drops. Our results highlight the importance of surface in controlling the motion of polyvalent BBRs.

Abrupt, Asynchronous Changes in Action Representations by Anterior Cingulate Cortex Neurons during Trial and Error Learning.

The ability to act on knowledge about the value of stimuli or actions factors into simple foraging behaviors as well as complex forms of decision-making. In striatal regions, action representations are thought to acquire value through a gradual (reinforcement-learning based) process. It is unclear whether this is also true for anterior cingulate cortex (ACC) where neuronal representations tend to change abruptly. We recorded from ensembles of ACC neurons as rats deduced which of 3 levers was rewarded each day. The rat's lever preferences changed gradually throughout the sessions as they eventually came to focus on the rewarded lever. Most individual neurons changed their responses to both rewarded and nonrewarded lever presses abruptly (<2 trials). These transitions occurred asynchronously across the population but peaked near the point where the rats began to focus on the rewarded lever. Because the individual transitions were asynchronous, the overall change at the population level appeared gradual. Abrupt transitions in action representations of ACC neurons may be part of a mechanism that alters choice strategies as new information is acquired.

DNA segregation under Par protein control.

The spatial organization of DNA is mediated by the Par protein system in some bacteria. ParB binds specifically to the parS sequence on DNA and orchestrates its motion by interacting with ParA bound to the nucleoid. In the case of plasmids, a single ParB bound plasmid is observed to execute oscillations between cell poles while multiple plasmids eventually settle at equal distances from each other along the cell's length. While the potential mechanism underlying the ParA-ParB interaction has been discussed, it remains unclear whether ParB-complex oscillations are stable limit cycles or merely decaying transients to a fixed point. How are dynamics affected by substrate length and the number of complexes? We present a deterministic model for ParA-ParB driven DNA segregation where the transition between stable arrangements and oscillatory behaviour depends only on five parameters: ParB-complex number, substrate length, ParA concentration, ParA hydrolysis rate and the ratio of the lengthscale over which the ParB complex stimulates ParA hydrolysis to the lengthscale over which ParA interacts with the ParB complex. When the system is buffered and the ParA rebinding rate is constant we find that ParB-complex dynamics is independent of substrate length and complex number above a minimum system size. Conversely, when ParA resources are limited, we find that changing substrate length and increasing complex number leads to counteracting mechanisms that can both generate or subdue oscillatory dynamics. We argue that cells may be poised near a critical level of ParA so that they can transition from oscillatory to fixed point dynamics as the cell cycle progresses so that they can both measure their size and faithfully partition their genetic material. Lastly, we show that by modifying the availability of ParA or depletion zone size, we can capture some of the observed differences in ParB-complex positioning between replicating chromosomes in B. subtilis cells and low-copy plasmids in E. coli cells.

Children's biobehavioral reactivity to challenge predicts DNA methylation in adolescence and emerging adulthood.

A growing body of research has documented associations between adverse childhood environments and DNA methylation, highlighting epigenetic processes as potential mechanisms through which early external contexts influence health across the life course. The present study tested a complementary hypothesis: indicators of children's early internal, biological, and behavioral responses to stressful challenges may also be linked to stable patterns of DNA methylation later in life. Children's autonomic nervous system reactivity, temperament, and mental health symptoms were prospectively assessed from infancy through early childhood, and principal components analysis (PCA) was applied to derive composites of biological and behavioral reactivity. Buccal epithelial cells were collected from participants at 15 and 18 years of age. Findings revealed an association between early life biobehavioral inhibition/disinhibition and DNA methylation across many genes. Notably, reactive, inhibited children were found to have decreased DNA methylation of the DLX5 and IGF2 genes at both time points, as compared to non-reactive, disinhibited children. Results of the present study are provisional but suggest that the gene's profile of DNA methylation may constitute a biomarker of normative or potentially pathological differences in reactivity. Overall, findings provide a foundation for future research to explore relations among epigenetic processes and differences in both individual-level biobehavioral risk and qualities of the early, external childhood environment.

Dense neural networks for predicting chromatin conformation.

BACKGROUND: DNA inside eukaryotic cells wraps around histones to form the 11nm chromatin fiber that can further fold into higher-order DNA loops, which may depend on the binding of architectural factors. Predicting how the DNA will fold given a distribution of bound factors, here viewed as a type of sequence, is currently an unsolved problem and several heterogeneous polymer models have shown that many features of the measured structure can be reproduced from simulations. However a model that determines the optimal connection between sequence and structure and that can rapidly assess the effects of varying either one is still lacking. RESULTS: Here we train a dense neural network to solve for the local folding of chromatin, connecting structure, represented as a contact map, to a sequence of bound chromatin factors. The network includes a convolutional filter that compresses the large number of bound chromatin factors into a single 1D sequence representation that is optimized for predicting structure. We also train a network to solve the inverse problem, namely given only structural information in the form of a contact map, predict the likely sequence of chromatin states that generated it. CONCLUSIONS: By carrying out sensitivity analysis on both networks, we are able to highlight the importance of chromatin contexts and neighborhoods for regulating long-range contacts, along with critical alterations that affect contact formation. Our analysis shows that the networks have learned physical insights that are informative and intuitive about this complex polymer problem.

A maximum-entropy model for predicting chromatin contacts.

The packaging of DNA inside a nucleus shows complex structure stabilized by a host of DNA-bound factors. Both the distribution of these factors and the contacts between different genomic locations of the DNA can now be measured on a genome-wide scale. This has advanced the development of models aimed at predicting the conformation of DNA given only the locations of bound factors-the chromatin folding problem. Here we present a maximum-entropy model that is able to predict a contact map representation of structure given a sequence of bound factors. Non-local effects due to the sequence neighborhood around contacting sites are found to be important for making accurate predictions. Lastly, we show that the model can be used to infer a sequence of bound factors given only a measurement of structure. This opens up the possibility for efficiently predicting sequence regions that may play a role in generating cell-type specific structural differences.

Operational Principles for the Dynamics of the In Vitro ParA-ParB System.

In many bacteria the ParA-ParB protein system is responsible for actively segregating DNA during replication. ParB proteins move by interacting with DNA bound ParA-ATP, stimulating their unbinding by catalyzing hydrolysis, that leads to rectified motion due to the creation of a wake of depleted ParA. Recent in vitro experiments have shown that a ParB covered magnetic bead can move with constant speed over a DNA covered substrate that is bound by ParA. It has been suggested that the formation of a gradient in ParA leads to diffusion-ratchet like motion of the ParB bead but how it forms and generates a force is still a matter of exploration. Here we develop a deterministic model for the in vitro ParA-ParB system and show that a ParA gradient can spontaneously form due to any amount of initial spatial noise in bound ParA. The speed of the bead is independent of this noise but depends on the ratio of the range of ParA-ParB force on the bead to that of removal of surface bound ParA by ParB. We find that at a particular ratio the speed attains a maximal value. We also consider ParA rebinding (including cooperativity) and ParA surface diffusion independently as mechanisms for ParA recovery on the surface. Depending on whether the DNA covered surface is undersaturated or saturated with ParA, we find that the bead can accelerate persistently or potentially stall. Our model highlights key requirements of the ParA-ParB driving force that are necessary for directed motion in the in vitro system that may provide insight into the in vivo dynamics of the ParA-ParB system.

Probing long-range interactions by extracting free energies from genome-wide chromosome conformation capture data.

BACKGROUND: A variety of DNA binding proteins are involved in regulating and shaping the packing of chromatin. They aid the formation of loops in the DNA that function to isolate different structural domains. A recent experimental technique, Hi-C, provides a method for determining the frequency of such looping between all distant parts of the genome. Given that the binding locations of many chromatin associated proteins have also been measured, it has been possible to make estimates for their influence on the long-range interactions as measured by Hi-C. However, a challenge in this analysis is the predominance of non-specific contacts that mask out the specific interactions of interest. RESULTS: We show that transforming the Hi-C contact frequencies into free energies gives a natural method for separating out the distance dependent non-specific interactions. In particular we apply Principal Component Analysis (PCA) to the transformed free energy matrix to identify the dominant modes of interaction. PCA identifies systematic effects as well as high frequency spatial noise in the Hi-C data which can be filtered out. Thus it can be used as a data driven approach for normalizing Hi-C data. We assess this PCA based normalization approach, along with several other normalization schemes, by fitting the transformed Hi-C data using a pairwise interaction model that takes as input the known locations of bound chromatin factors. The result of fitting is a set of predictions for the coupling energies between the various chromatin factors and their effect on the energetics of looping. We show that the quality of the fit can be used as a means to determine how much PCA filtering should be applied to the Hi-C data. CONCLUSIONS: We find that the different normalizations of the Hi-C data vary in the quality of fit to the pairwise interaction model. PCA filtering can improve the fit, and the predicted coupling energies lead to biologically meaningful insights for how various chromatin bound factors influence the stability of DNA loops in chromatin.

Concordant and discordant DNA methylation signatures of aging in human blood and brain.

BACKGROUND: DNA methylation is an epigenetic mark that balances plasticity with stability. While DNA methylation exhibits tissue specificity, it can also vary with age and potentially environmental exposures. In studies of DNA methylation, samples from specific tissues, especially brain, are frequently limited and so surrogate tissues are often used. As yet, we do not fully understand how DNA methylation profiles of these surrogate tissues relate to the profiles of the central tissue of interest. RESULTS: We have adapted principal component analysis to analyze data from the Illumina 450K Human Methylation array using a set of 17 individuals with 3 brain regions and whole blood. All of the top five principal components in our analysis were associated with a variable of interest: principal component 1 (PC1) differentiated brain from blood, PCs 2 and 3 were representative of tissue composition within brain and blood, respectively, and PCs 4 and 5 were associated with age of the individual (PC4 in brain and PC5 in both brain and blood). We validated our age-related PCs in four independent sample sets, including additional brain and blood samples and liver and buccal cells. Gene ontology analysis of all five PCs showed enrichment for processes that inform on the functions of each PC. CONCLUSIONS: Principal component analysis (PCA) allows simultaneous and independent analysis of tissue composition and other phenotypes of interest. We discovered an epigenetic signature of age that is not associated with cell type composition and required no correction for cellular heterogeneity.

Bacterial motion in narrow capillaries.

Motile bacteria often have to pass through small tortuous pores in soil or tissue of higher organisms. However, their motion in this prevalent type of niche is not fully understood. Here, we modeled it with narrow glass capillaries and identified a critical radius (Rc) for bacterial motion. Near the surface of capillaries narrower than that, the swimming trajectories are helices. In larger capillaries, they swim in distorted circles. Under non-slip condition, the peritrichous Escherichia coli swam in left-handed helices with an Rc of ~10 µm near glass surface. However, slipping could occur in the fast monotrichous Pseudomonas fluorescens, when a speed threshold was exceeded, and thus both left-handed and right-handed helices were executed in glass capillaries. In the natural non-cylindrical pores, the near-surface trajectories would be spirals and twisted loops. Engaging in such motions reduces the bacterial migration rate. With a given pore size, the run length and the tumbling angle of the bacterium determine the probability and duration of their near-surface motion. Shear flow and chemotaxis potentially enhance it. Based on this observation, the puzzling previous observations on bacterial migration in porous environments can be interpreted.

Confinement-dependent localization of diffusing aggregates in cellular geometries.

Confinement has a strong influence on diffusing nano-sized clusters. In particular, biomolecular aggregates within the shell-like confining space of a bacterial cell have been shown to display a variety of localization patterns, from being midcell to the poles. How does the confining space determine where the aggregate will localize? Here, using Monte Carlo simulations we have calculated the equilibrium spatial distribution of fixed-sized clusters diffusing in spherocylindrical shells. We find that localization to the poles depends strongly on shell thickness and the size of the cluster. Compared to being at midcell, polar clusters can be more bent and hence have higher energy, but they also can have a greater number of defects and hence have more entropy. Under certain conditions this can lead to polar clusters having a lower free energy than at midcell, favoring localization to the poles. Our findings suggest possible localization selection mechanisms within shell-like geometries that can arise purely from cluster confinement.

Localization of aggregating proteins in bacteria depends on the rate of addition.

Many proteins are observed to localize to specific subcellular regions within bacteria. Recent experiments have shown that proteins that have self-interactions that lead them to aggregate tend to localize to the poles. Theoretical modeling of the localization of aggregating protein within bacterial cell geometries shows that aggregates can spontaneously localize to the pole due to nucleoid occlusion. The resulting polar localization, whether it be to a single pole or to both was shown to depend on the rate of protein addition. Motivated by these predictions we selected a set of genes from Escherichia coli, whose protein products have been reported to localize when tagged with green fluorescent protein (GFP), and explored the dynamics of their localization. We induced protein expression from each gene at different rates and found that in all cases unipolar patterning is favored at low rates of expression whereas bipolar is favored at higher rates of expression. Our findings are consistent with the predictions of the model, suggesting that localization may be due to aggregation plus nucleoid occlusion. When we expressed GFP by itself under the same conditions, no localization was observed. These experiments highlight the potential importance of protein aggregation, nucleoid occlusion and rate of protein expression in driving polar localization of functional proteins in bacteria.

Insulators recruit histone methyltransferase dMes4 to regulate chromatin of flanking genes.

Chromosomal domains in Drosophila are marked by the insulator-binding proteins (IBPs) dCTCF/Beaf32 and cofactors that participate in regulating long-range interactions. Chromosomal borders are further enriched in specific histone modifications, yet the role of histone modifiers and nucleosome dynamics in this context remains largely unknown. Here, we show that IBP depletion impairs nucleosome dynamics specifically at the promoters and coding sequence of genes flanked by IBP binding sites. Biochemical purification identifies the H3K36 histone methyltransferase NSD/dMes-4 as a novel IBP cofactor, which specifically co-regulates the chromatin accessibility of hundreds of genes flanked by dCTCF/Beaf32. NSD/dMes-4 presets chromatin before the recruitment of transcriptional activators including DREF that triggers Set2/Hypb-dependent H3K36 trimethylation, nucleosome positioning, and RNA splicing. Our results unveil a model for how IBPs regulate nucleosome dynamics and gene expression through NSD/dMes-4, which may regulate H3K27me3 spreading. Our data uncover how IBPs dynamically regulate chromatin organization depending on distinct cofactors.

Chromatin immunoprecipitation indirect peaks highlight long-range interactions of insulator proteins and Pol II pausing.

Eukaryotic chromosomes are partitioned into topologically associating domains (TADs) that are demarcated by distinct insulator-binding proteins (IBPs) in Drosophila. Whether IBPs regulate specific long-range contacts and how this may impact gene expression remains unclear. Here we identify "indirect peaks" of multiple IBPs that represent their distant sites of interactions through long-range contacts. Indirect peaks depend on protein-protein interactions among multiple IBPs and their common cofactors, including CP190, as confirmed by high-resolution analyses of long-range contacts. Mutant IBPs unable to interact with CP190 impair long-range contacts as well as the expression of hundreds of distant genes that are specifically flanked by indirect peaks. Regulation of distant genes strongly correlates with RNAPII pausing, highlighting how this key transcriptional stage may trap insulator-based long-range interactions. Our data illustrate how indirect peaks may decipher gene regulatory networks through specific long-range interactions.

Distinct DNA methylation patterns of cognitive impairment and trisomy 21 in Down syndrome.

BACKGROUND: The presence of an extra whole or part of chromosome 21 in people with Down syndrome (DS) is associated with multiple neurological changes, including pathological aging that often meets the criteria for Alzheimer's Disease (AD). In addition, trisomies have been shown to disrupt normal epigenetic marks across the genome, perhaps in response to changes in gene dosage. We hypothesized that trisomy 21 would result in global epigenetic changes across all participants, and that DS patients with cognitive impairment would show an additional epigenetic signature. METHODS: We therefore examined whole-genome DNA methylation in buccal epithelial cells of 10 adults with DS and 10 controls to determine whether patterns of DNA methylation were correlated with DS and/or cognitive impairment. In addition we examined DNA methylation at the APP gene itself, to see whether there were changes in DNA methylation in this population. Using the Illumina Infinium 450 K Human Methylation Array, we examined more than 485,000 CpG sites distributed across the genome in buccal epithelial cells. RESULTS: We found 3300 CpGs to be differentially methylated between the groups, including 495 CpGs that overlap with clusters of differentially methylated probes. In addition, we found 5 probes that were correlated with cognitive function including two probes in the TSC2 gene that has previously been associated with Alzheimer's disease pathology. We found no enrichment on chromosome 21 in either case, and targeted analysis of the APP gene revealed weak evidence for epigenetic impacts related to the AD phenotype. CONCLUSIONS: Overall, our results indicated that both Trisomy 21 and cognitive impairment were associated with distinct patterns of DNA methylation.

A model for cell population size control using asymmetric division.

In multicellular organisms one can find examples where a growing tissue divides up until some final fixed cell number. Asymmetric division plays a prevalent feature in tissue differentiation in these organisms, where the daughters of each asymmetric division inherit unequal amounts of a fate determining molecule and as a result follow different developmental fates. In some tissues the accumulation or decrease of cell cycle regulators acts as an intrinsic timing mechanism governing proliferation. Here we present a minimal model based on asymmetric division and dilution of a cell-cycle regulator that can generate any final population size that might be needed. We show that within the model there are a variety of growth mechanisms from linear to non-linear that can lead to the same final cell count. Interestingly, when we include noise at division we find that there are special final cell population sizes that can be generated with high confidence that are flanked by population sizes that are less robust to division noise. When we include further perturbations in the division process we find that these special populations can remain relatively stable and in some cases even improve in their fidelity.

Abrupt changes in the patterns and complexity of anterior cingulate cortex activity when food is introduced into an environment.

Foraging typically involves two distinct phases, an exploration phase where an organism explores its local environment in search of needed resources and an exploitation phase where a discovered resource is consumed. The behavior and cognitive requirements of exploration and exploitation are quite different and yet organisms can quickly and efficiently switch between them many times during a foraging bout. The present study investigated neural activity state dynamics in the anterior cingulate sub-region of the rat medial prefrontal cortex (mPFC) when a reliable food source was introduced into an environment. Distinct and largely independent states were detected using a Hidden Markov Model (HMM) when food was present or absent in the environment. Measures of neural entropy or complexity decreased when rats went from exploring the environment to exploiting a reliable food source. Exploration in the absence of food was associated with many weak activity states, while bouts of food consumption were characterized by fewer stronger states. Widespread activity state changes in the mPFC may help to inform foraging decisions and focus behavior on what is currently most prominent or valuable in the environment.